Statistics in Medicine

Estimating the correlation between subject-specific rates of change is a common objective in longitudinal studies. A natural approach is to fit ordinary least squares regressions per subject and correlate the resulting slope estimates directly. Under complete data, this slopes-first estimator is algebraically equivalent to a contrast-based plug-in estimator derived from the joint variance-covariance matrix of the two signals across time. Under independent dropout, however, the slopes-first approach implicitly applies heterogeneous weights across missing-data patterns, introducing bias whenever the cross-covariance between the two signals is not constant over time. We formalise this result and propose an unbiased alternative that estimates the joint covariance matrix before applying the contrast. We further extend our results to include covariates as well. Simulations confirm that the slopes-first estimator is biased under non-constant cross-covariance while the proposed alternative remains unbiased. Further, we show in a real-world example that the two estimators can lead to different scientific conclusions and recommend that researchers working with longitudinal data inspect the missing data pattern to make an informed decision about the estimator of choice.

Background: Cost-effectiveness analyses (CEA) built on simulation models are usually delivered through project-specific scripts. Each new model, and often each new question, rebuilds the same scaffolding: parameter management, calibration, base-case analysis, sensitivity analysis and export. This is costly, and is where traceability between an estimate and the model version producing it is lost. We present THALASSA, a web application providing this infrastructure once, for any model satisfying a small programming contract. Methods: A simulation model is connected by one file implementing a fixed contract of R functions that determine its behaviour. THALASSA exposes the complete CEA workflow in a browser: a generated overview of the connected model, parameter editing, calibration, base-case comparison with the efficiency curve, and deterministic and probabilistic sensitivity analysis at a willingness-to-pay threshold, each with timestamped file export carrying the parameter values it was run with. Every session is bound to a repository version, so results carry their generating model version. Results: We demonstrate THALASSA on two public models of contrasting difficulty: a simple tutorial model and a generic cancer screening model with two competing natural-history pathways and a 21-dimensional calibration deliberately built to be more complex. Both ran on the same unmodified setup, implementing the same contract, with no code changes. Conclusions: Separating analysis infrastructure from model logic makes the CEA workflow reusable across projects and reproducible by construction. A modelling team can move from a versioned repository to audit-ready, exportable CEA outputs without writing analysis code. The integration effort does not grow with model size and one installation serves unrelated projects.

Background and Objective: Adding belzutifan to adjuvant pembrolizumab improved disease-free survival (DFS) in LITESPARK-022, but its durability beyond the 1-year treatment is unknown. Methods: Following a registered plan, we reconstructed individual patient data from published Kaplan-Meier curves of ten phase 3 adjuvant renal cell carcinoma trials. Within each trial, we estimated restricted mean survival time differences ({Delta}RMST), hazard ratios (HRs) before and after planned treatment ended, and DFS differences. Key Findings and Limitations: In LITESPARK-022 (central review), the 24-month {Delta}RMST was 1.34 months (95% confidence interval 0.73-1.95), and the belzutifan arm remained ahead after month 12 ({Delta}RMST increment 0.74 months/yr, 0.32-1.19). The HR was 0.55 (0.43-0.71) before month 12 and 1.16 (0.87-1.54) afterwards (post/pre ratio 2.10); investigator assessment showed similar attenuation. Similar ratios occurred in KEYNOTE-564 at its first interim analysis (2.02), whose benefit persisted, and in SORCE and PROTECT (2.09, 2.59), whose benefit faded. From months 12 to 30 (post hoc), the DFS difference changed by -3.7 percentage points (-8.1 to 0.8) with central review and 0.6 (-4.0 to 5.0) with investigator assessment. Data were reconstructed, and follow-up was short. Conclusions and Clinical Implications: Relative hazard reductions with belzutifan were greater before month 12; durability of the absolute DFS advantage is uncertain and assessment-dependent, so single-HR projections should be avoided.

People experiencing homelessness (PEH) have substantially lower quality-adjusted life expectancy and substantially higher healthcare spending compared to housed populations. California has the largest population of PEH in the United States, and an estimated ten percent of PEH have co-occurring opioid use disorder (OUD). We developed a lifetime, individual-level microsimulation of the 18,928 PEH with OUD in California to evaluate the health and economic effects of expanding access to housing through a Housing First approach, which does not require OUD treatment or abstinence. Increases in permanent housing and shelter access were based on Los Angeles County's Measure A goals, and we accounted for both housing and healthcare costs. Expanding housing access increased mean age at death by 0.34 years, reduced fatal overdoses for the cohort by 39, and produced an incremental cost-effectiveness ratio of $31,083 per quality-adjusted life-year gained, well below conventional willingness-to-pay thresholds. We provide evidence that housing is a high-value health intervention for PEH with OUD in California. Policymakers should support Housing First programs as part of integrated responses to the intertwined crises of homelessness and opioid use.

The efficacy of current antiretroviral interventions for prevention of vertical HIV transmission is both a triumph and a challenge for future efficacy studies of novel agents, including long-acting antiretroviral agents. Given the effectiveness of current antiretroviral interventions, standard randomized controlled trials to assess novel interventions targeted towards further reducing vertical HIV transmission may be infeasible due to required sample sizes. To address this study design challenge, we propose a novel augmented clinical trial design, where a randomized trial with an active comparator is augmented with a concurrent observational cohort to supplement the standard-of-care arm. The proposed design and analysis can overcome sampling biases, and an associated diagnostic can evaluate potential biases from unmeasured covariates. The proposed trial design is illustrated with an example and studied further using a simulation study. In the illustrative example, an augmented trial with a 3:1 allocation ratio in the trial component yielded an unbiased treatment effect estimate, as compared to the trial component alone, and the standard error was reduced substantially. Although the total sample size of the augmented trial was larger than a similarly powered 1:1 standard trial, the sample size of the augmented trial's trial component was half the size of the standard trial and within the range of historically similar trials, thereby providing both feasibility and potential cost advantages. Although limitations with this design remain, many of these limitations can be mitigated by careful study planning and coordination. Augmented trials provide a potentially applicable solution to studying the efficacy of long-acting antiretrovirals for infant postnatal prophylaxis against HIV.

The aim of this study was to develop a framework for assessing the incremental value of decisions on drug amount, timing, input structure, cessation and restart. The LOGOS (Longitudinal Optimization of Global Outcomes and Strategies) approach combines matched nested regimen classes, mechanism-specific benchmarks and evidence requirements in a single longitudinal comparison. Thus, the contribution of additional regimen freedom is assessed separately from improvements due to different resources, constraints or evaluation horizons. Nested feasible classes have nondecreasing objective suprema, whereas a strict gain requires a feasible witness above a valid bound for the smaller class. The approach is illustrated through analytical examples of periodic smoothing, matched desensitization, relaxed exposure allocation, finite-course loading and safety after cessation. In a normalized finite Emax model, integrated benefit increased from 0.629606 for three equal boluses to 0.672945 after optimization of amounts, 0.677043 for an amount-and-time candidate and 0.679553 for a globally certified loading-maintenance-coast solution. The adjacent observed gains were 6.88%, 0.61% and 0.37%, respectively, and all schedules satisfied the full-tail injury bound. In addition, a periodic factorial example separated response curvature from a joint recovery and repair process, including a null case and results dependent on feasibility constraints. The LOGOS framework allows assessment of whether additional regimen flexibility improves benefit, provides no added value or restores feasibility. These dimensionless examples support conditional mathematical and model-based conclusions; clinical superiority and safe restart in patients have not been established.

Polygenic scores (PGSs) have progressed toward clinical utility in recent years, demonstrating odds ratios that are comparable with both lifestyle-related and monogenic risk factors for several complex diseases. However, PGS implementation relies on developing analytically and clinically valid assays. Methods for developing polygenic models (PGMs) that incorporate tens of thousands to millions of variants require either genome-wide chip genotyping combined with imputation or whole genome sequencing, limiting their portability. Additionally, variants may be in regions where the local genomic sequence makes it difficult to assay, and multi-allelic or strand ambiguous variants can make bioinformatics more complicated. Because of these challenges, attempts to develop a reliable assay for published PGMs can produce multiple variations from the same underlying model, hindering routine clinical deployment and incorporation into multi-factorial risk prediction models. In this analysis, we present an adjustment to the previously published select and shrink with summary statistics (S4) PGM method designed to optimize performance for realistic clinical use by restricting to a pre-specified number of available SNPs that are well-captured by genotyping arrays. We compare the performance of this oligogenic S4 method (S4-Oligo) to currently published PGSs as well as the full S4 model for prediction of breast, colorectal, ovarian, and prostate cancers in the All of Us Research Program, Biobank Japan, and UK Biobank. We find that S4-Oligo matches or exceeds the accuracy of published PGMs at comparable model sizes while retaining between 88% and 99% of the performance of the original S4 model when limited to 1,000 variants.

Personalized medicine aims to tailor treatments to individual patients, but its implementation is limited by scarce clinical resources and heterogeneous treatment effects. Electronic health records now offer longitudinal data that can support data-driven therapy selection. We focus on SGLT2 inhibitors, medications with broad metabolic and cardiovascular benefits, whose clinical impact varies significantly based on individual patient profiles. We propose a counterfactual Structural Causal Model (SCM) to estimate individual treatment effects (ITE) of SGLT2 inhibitors from observational data. The model assumes counterfactual stability and uses a Gumbel-Max sampling mechanism to capture stochastic treatment assignment while enabling differentiable inference. Optimal treatments are defined by maximizing the counterfactual value of future diagnostic outcomes. We evaluate the Counterfactual Gumbel-Max SCM using a doubly robust value estimator and show that, across diagnostic categories, it outperforms real-world prescriptions in predicted treatment value, highlighting subpopulations where SGLT2 inhibitors yield superior outcomes. By integrating causal inference with electronic health data, this framework provides individualized counterfactual reasoning about SGLT2 treatment effects and can support clinicians in identifying optimal therapies for type 2 diabetes and comorbid conditions.

Importance: Coccidioidomycosis, also known as Valley fever, is a serious fungal disease caused by inhaling Coccidioides spores aerosolized from the soil. Despite rising public health concern, the distribution and incidence of coccidioidomycosis in the US remain poorly characterized as not all states are required to report cases, and case under-ascertainment is presumed to be high, particularly outside of the endemic region. Objective: To estimate monthly, county-level coccidioidomycosis incidence across the US while accounting for widespread under-ascertainment. Design, Setting, and Participants: We integrated multisource data, including human cases diagnosed in electronic health records, serological testing of dogs, and ecological and sociodemographic predictors, in a machine learning model estimating monthly, county-level coccidioidomycosis incidence. The model was trained on surveillance data from Arizona and California from 2013 to 2023, then applied to predict incidence in states across the contiguous US over the same time period. Main Outcomes and Measures: Incidence of coccidioidomycosis at the monthly county level. Results: The model predicted coccidioidomycosis incidence from 2013 to 2023 for 3,108 counties from all states in the contiguous US. The model identified 60 counties in non-reporting states where predicted incidence was [≥]2 cases per 100,000 population (considered moderate to high incidence). In states with mandatory reporting outside of Arizona and California, the model predicted county-level incidence rates that were 28% higher (95% bootstrap interval, 0.1%-69.8%) than those captured in national surveillance, equating to a total of 608 potentially missed cases per year (95% bootstrap interval, 340-964). Conclusions and Relevance: In this cross-sectional study, the model quantified case under-ascertainment for coccidioidomycosis relative to Arizona and California, and highlighted the potential for alternative data sources to fill in gaps in mandatory reporting. Estimating incidence while accounting for under-ascertainment and non-reporting may help efforts to improve public and physician awareness and to detect emergence in new regions.

We develop an age- and immunity-structured COVID-19 model with setting-specific contacts, calibrated to Ontario incidence and hospital admissions before Omicron. We simulate counterfactual vaccination strategies and contact restrictions in schools, workplaces, and other community settings, evaluating projected severe-case demand against an illustrative healthcare-capacity benchmark. Vaccination combined with workplace and community contact restrictions keeps simulated demand within capacity. In contrast, relaxing restrictions in either setting produces large surges and repeated capacity threats, even when school restrictions remain. School-only changes produce smaller effects on severe disease. These findings highlight the value of vaccination and setting-specific contact restrictions that prioritize workplace and community settings for pathogens with adult-concentrated severe outcomes.

Background: TB vaccine development is hampered by the natural history of M. tuberculosis: most infected people rarely progress to disease, so trials must be very large and long to accrue sufficient disease endpoints. Current trials enrol IGRA-positive participants (IGRA: Interferon-Gamma Release Assays), but IGRAs remain positive for many years and cannot distinguish recent from remote infection. Tests for recent infection are now in development, and we argue they could make trials substantially more efficient. Here, we investigate how screening participants prior to enrolment with a test for recent M.tb infection could impact the total number of enrollees required for an adequately-powered TB vaccine trial. Methods: Using a statistical model of vaccine trial recruitment, we estimated the number needed to contact, test, and enrol to achieve an 'effective' trial (defined as 50 disease endpoints), under IGRA-gated enrolment (life-long duration of positivity) vs. enrolment with a test for recent infection. We considered five recruitment strategies: non-test-based enrolment for ages 15-49, ages 12-34, and household contacts of recent cases; and test-based enrolment using IGRA, and a test for recent infection with two-year duration of positivity. The model explored trial recruitment in South Africa, Kenya, and India. Results: Compared with an extended-positivity test, a test for recent infection reduced the median number needed to enrol by 4.2- to 6-fold under favourable conditions. It also outperformed non-test-based strategies, including enrolling household contacts. This comes at the cost of substantially more testing, but crude estimates suggest recency-based enrolment would still be cheaper overall provided the new test does not cost more than roughly 5-10 times an IGRA. Conclusions: Diagnostics for recent M.tb infection have the potential to improve vaccine trial efficiency, but more research is now needed to understand test characteristics and to evaluate their use in this setting. Our Viewpoint sets out the test attributes that would need to be established, such as time to positivity, duration of positivity, sensitivity, and specificity, before such tests could be used in vaccine trials.

Background. Fluid collections are common after kidney transplantation, but their prevalence at scheduled imaging and their management are inconsistent from canter to canter. We describe the cross sectional prevalence and management of perigraft collections, including an exploratory deepdrain subgroup. Methods. In this retrospective single center cohort, 146 consecutive adult kidney transplant recipients (Cleveland Clinic Weston Hospital, July 2024 to June 2025) were imaged on postoperative day 0/1 and at approximately 1 and 3 months. Routine deep drainage was avoided (n=124), 22 clinically selected recipients with perioperative anticoagulation received a deep drain. The primary outcome was cross sectional prevalence of a perigraft collection at each visit. Secondary outcomes were management among patients with a collection. Drain group comparisons were exploratory. Results. Collections were present in 29 of 146 recipients at postoperative day 0/1 (19.9%; 95% CI, 14.2 to 27.1), 71 of 146 at 1 month (48.6%; 95% CI, 40.7 to 56.7), and 34 of 146 at 3 months (23.3%; 95% CI, 17.2 to 30.8). Prevalence in the no drain group was 21.0%, 50.0%, and 23.4%, respectively; corresponding prevalence in the deep drain group was 13.6%, 40.9%, and 22.7%. Between group prevalence comparisons were not statistically significant at any time point. BMI was higher in the deep-drain group than in the no-drain group (30.04 vs 27.1 kg/m2; P=0.01). Among patients with a 1-month collection, management did not differ significantly between groups (any invasive management, 11.1% vs 16.1%; P=1.00). Among patients with a 3month collection, invasive management was more common in the deep drain group than in the no drain group (80.0% vs 13.8%; P=0.007). This comparison included only five deep-drain recipients and was confounded by clinical selection, perioperative anticoagulation, and BMI imbalance. Conclusions. Perigraft collection prevalence was highest at 1 month and lower at 3 months in this predominantly no drain cohort. Most collections were managed conservatively, and invasive treatment was uncommon in the no-drain group. The deep drain subgroup provides clinically relevant information about high risk recipients selected for drainage but cannot support causal conclusions regarding drain benefit or harm. Patient level linkage is required to distinguish persistence, resolution, and new collection formation.

Retention in antiretroviral therapy (ART) programmes remains a major public health challenge in HIV care, with loss to follow-up (LTFU) and death threatening treatment success. Because these outcomes are competing events, conventional survival methods may yield biased estimates of their incidence and associated risk factors. To estimate the cumulative incidence of LTFU and mortality and identify factors associated with these outcomes among adults receiving ART in Northern Ghana using a competing-risks framework. We conducted a retrospective cohort study of 329 adults living with HIV who initiated ART at Nalerigu Baptist Medical Centre, Ghana, between January 2013 and December 2019. Fine--Gray subdistribution hazard models were used to evaluate predictors of LTFU and death while accounting for competing risks. Cause-specific Cox proportional hazards models were fitted as a sensitivity analysis. The incidence rates were 2.2 per 100 person-months (95% CI: 1.9--2.6) for LTFU and 0.74 per 100 person-months (95% CI: 0.57--0.98) for death. At 36 months, the cumulative incidence was 50.2% for LTFU and 16.0% for death. Patients residing closer to the clinic had a lower cumulative incidence of LTFU than those living farther away (SHR = 0.68, 95% CI: 0.49--0.93), whereas higher baseline weight was associated with increased LTFU (SHR = 1.02 per kg, 95% CI: 1.00--1.04). Mortality was lower among patients receiving efavirenz-based regimens (SHR = 0.27, 95% CI: 0.15--0.50), females (SHR = 0.39, 95% CI: 0.18--0.85), and those with higher baseline weight (SHR = 0.93 per kg, 95% CI: 0.90--0.96). Findings were broadly consistent in cause-specific Cox models. LTFU was substantially more frequent than mortality among adults receiving ART in Northern Ghana. Travel distance to the treatment facility was an important barrier to retention, while ART regimen, sex, and baseline weight were associated with mortality. Expanding access to ART services through decentralized and community-based models may improve retention and support efforts to achieve HIV control targets.

Background: Noninvasive biomarkers capable of detecting kidney allograft injury and distinguishing rejection from infection remain an important clinical need. Torque teno virus (TTV) reflects host immune status, but its utility in urine has not been well established. We evaluated urinary TTV, total cell-free DNA (cfDNA), cfDNA fragmentation, and BK polyomavirus (BKV) as potential biomarkers of kidney transplant rejection and BK polyomavirus-associated nephropathy (BKPyVAN). Methods: Urine samples from 334 kidney transplant recipients were analyzed in association with kidney biopsies performed within 30 days of sample collection. TTV and cfDNA were measured by multiplex quantitative PCR (qPCR), with a subset analyzed by digital PCR (dPCR), and BKV was quantified by qPCR. Associations with biopsy results were evaluated using group comparisons and nominal logistic regression. Changes associated with outcomes of consecutive biopsies were assessed by matched pairs analysis. Results: Among 544 unique urine samples, the urinary TTV/cfDNA ratio was significantly lower in rejection than in non-rejection samples (P<0.0001). Logistic regression adjusted for demographic factors confirmed an inverse association between the qPCR-derived TTV/cfDNA metric and rejection (adjusted odds ratio per delta Cq, 0.91; 95% CI, 0.88-0.94; P<0.00001). Significant differences were observed for all types of rejection. Association was confirmed by dPCR: higher TTV was associated with lower rejection odds, whereas higher cfDNA was associated with greater rejection odds. A nominal logistic model incorporating dPCR measurements achieved an AUROC of 0.76 (95% CI, 0.69-0.82). Matched pairs analysis showed a significant decrease in TTV/cfDNA when rejection developed and an increase following recovery from rejection, with no significant changes when biopsy status remained unchanged. BKV levels were strongly associated with BKPyVAN. No correlation between TTV and BKV was observed. Nominal logistic regression integrating TTV, cfDNA, and BKV measurements yielded AUROCs of 0.75, 0.97, and 0.73 for rejection, BKPyVAN, and non-rejection, respectively. Conclusions: A reduced urinary TTV/cfDNA ratio is associated with biopsy-defined kidney transplant rejection, and changes in the ratio correlate with transitions in biopsy results. Combined TTV, cfDNA, and BKV measurements showed differential associations with rejection and BKPyVAN, supporting further evaluation as complementary urinary biomarkers.

Importance: Preventing substance use disorder (SUD) is a top public health priority, but clinical practice currently lacks validated prognostic tools to identify risk before the onset of heavy use/problems. Objective: To determine whether models can predict SUD using genetic, behavioral, and environmental factors and to evaluate performance across substance types (alcohol, cannabis, tobacco), sex, and genetic similarity groups. Design: This prognostic study developed and validated a multi-domain risk model using data from four longitudinal cohorts (The National Longitudinal Study of Adolescent to Adult Health, Avon Longitudinal Study of Parents and Children, The Collaborative Study on the Genetics of Alcoholism, and FinnTwin12) with follow-up from early adolescence into young adulthood. Setting: Population-based and family-based samples from the United States, United Kingdom, and Finland. Participants: The analysis included participants with available genomic and behavioral data. The cohorts represented a spectrum of risk, ranging from nationally representative community samples to high-risk populations. Exposure(s) (for observational studies): Predictors included polygenic scores for substance use-related phenotypes and 10 established early-life behavioral and environmental risk factors. Main Outcome(s) and Measure(s): The outcomes were a diagnosis of moderate-to-severe SUD based on the Diagnostic and Statistical Manual of Mental Disorders, Fifth Edition (DSM-5) criteria. Model performance was evaluated using the area under the receiver operating characteristic curve (AUC), sensitivity, specificity, and positive and negative predictive value. Results: A total of 16,362 participants were included in the study (54.7% female; mean [SD] age 26.73 [3.36]). Genetic similarity group composition was 70.4% European, 16.8% African, 9.8% admixed American, 2.2% East Asian, and 0.7% South Asian. SUD prevalence was 22.0%. The integrated model achieved an AUC of 0.78 for predicting any SUD, with a sensitivity of 69% and a positive predictive value of 45%. Performance was comparable for specific substances (AUC range 0.74-0.82) and across genetic similarity groups (AUC range 0.61-0.82). The model demonstrated higher AUC in females compared with males. Conclusions and Relevance: Findings highlight the potential utility of an SUD prediction model to inform clinicians and patients about risk for developing a SUD.

Background: Current Lynch syndrome (LS) prediction models do not account for detailed family structure and tumour mismatch repair (MMR) status. We developed LynchRisk, a Bayesian model that estimates individualized LS probability from pedigree, cancer history and tumour MMR test results. Methods: LynchRisk combines Mendelian inheritance, age- and sex-specific cancer hazards, affected and unaffected relatives, and detailed MMR status. Parameters were derived from the literature. Results: In two simulated cases of probands with LS-associated cancers, family structure and MMR status markedly changed LS probability estimates, with implications for germline testing decisions. LynchRisk also identified situations in which a negative germline result might be considered sufficient to rule out LS or, conversely, prompt further investigations such as somatic MMR or long-read sequencing. Discussion: LynchRisk gives individualized LS probability by combining clinical, family, and tumour data at different stages of the diagnostic process. Further work will focus on recalibrating cancer hazard functions, validating model performance through simulations, and developing an R package for use by interested research groups.

Purpose: Accurately capturing a patient's health journey is essential for ensuring the validity of scientific and public health studies. Electronic health records (EHR) and administrative claims are commonly used in healthcare research, albeit separately. Previous research has demonstrated that EHR and claims capture different aspects of patient history and care, potentially introducing bias when used alone. In this study, we leveraged the largest public US EHR resource, the National Clinical Cohort Collaborative (N3C), alongside Centers for Medicare & Medicaid Services (CMS) claims data to compare their information coverage and agreement. Methods: We integrated N3C EHR and CMS data to create one cohesive dataset and quantified the proportion of key variables present in EHR alone, CMS alone, and the overlap present in both data sources. Overlap in information coverage was quantified pre- and post-landmark diagnosis of COVID-19. To evaluate the impact of data coverage differences, we conducted a case study evaluating the effect of COVID-19 vaccination on ischemic stroke using both data sources alone and in combination. Results: CMS data helped identify many individuals that were duplicated among EHR identifiers. Examination of the overlap showed that CMS data more completely captured conditions, while EHR data more completely captured medications after COVID-19 diagnosis. Some variables exhibited marked differences over time after landmark COVID-19 diagnosis. Using target trial emulation and self-controlled case series methods, we found that data capture variations altered effect estimates: EHR data alone yielded lower risk of ischemic stroke after vaccination in the target trial emulation (odds ratio 0.47, 95% CI 0.28-0.79), with a smaller, non-significant reduction in the self-controlled case series (incidence rate ratio 0.81, 95% CI 0.64-1.02), whereas analyses using CMS data alone or both sources combined yielded estimates close to the null in both designs. Conclusions: Results indicated that when using either a national database combining multiple EHR sources, or CMS claims data alone, patients' medical histories are less completely captured compared to when these data sources are combined. Case study results highlighted the value of combining data sources for not only yielding fuller patient health histories and valid scientific results, but for identifying and explaining biases that may emerge in research studies that rely on a single data source.

This paper examines the impact of Cancer Care Ontario's regionalization of lung cancer surgery on blood transfusion use for thoracic surgery at The Ottawa Hospital. Using aggregate data scraped from the paper ''Changing trends in blood transfusion: an analysis of 244,013 hospitalizations'' by Shehata et al. (2014) we estimated the effect regionalization on the log mean number of blood transfusions among thoracic surgeons (who treat lung cancer), utilizing other surgical specialties as controls. We found that the average number of blood transfusions increased significantly among thoracic surgeons following the opening of a new cancer care centre at the hospital in December 2009, while control specialties showed no parallel change. Contextualizing these results in a modified ''population level'' treatment threshold model, it is argued that the main driver of this increase in blood transfusions among thoracic surgeons was due to the influx of less skilled thoracic surgeons to The Ottawa Hospital, and the closing of thoracic surgery units at smaller hospitals which was part of the broader regionalization process.

Background: Iron deficiency (ID) and iron deficiency anaemia (IDA) have global and severe health implications, and disproportionately impact women of reproductive age. Current methods of supplementing diets with iron are challenged by poor adherence and palatability, reducing effectiveness of these treatment strategies. The aim of this study is to evaluate the efficacy of a novel home beverage fortification with a mineral-enriched powder, designed to enhance bioavailability and ease of use, for the treatment of ID in women of reproductive age. Methods: This single-centre, randomised, double-blinded placebo-controlled trial will enrol 150 eligible biologically female participants with ID. Participants will be randomly assigned to the active intervention (iron/zinc-enriched powder) or control group (placebo powder) over 6 months. Primary outcomes will include change in serum ferritin (SF) concentrations and proportion of participants who remain iron deficient at 6-months. Secondary outcomes will include change in SF and haemoglobin (Hb) concentrations across timepoints, and levels of metabolites and key circulating pro-inflammatory biomarkers, changes in dietary intakes, self-reported perceptions of health, barriers to use of the powder, adherence and palatability, safety outcomes, and cost-effectiveness. Exploratory outcomes will assess the gut microbial metagenome and metabolome. Discussion: Fortification, an effective approach for the prevention of ID/IDA, has advantages over some methods of supplementation because it uses low levels of iron, minimising unpleasant side effects, does not risk toxicity, and has high acceptability and compliance rates. This mineral-enriched powder uses electrolytic iron and zinc sulphate monohydrate of a small particle size, which is less irritating to the gut, and facilitates bioavailability. Findings from the study will be shared with those affected by and involved in treating these conditions, as a novel intervention that is both effective and acceptable among women of reproductive age would help to improve treatment compliance, and ultimately, reduce ID/IDA burden. Trial registration: This trial is registered at clinicaltrials.gov (NCT05990166): https://clinicaltrials.gov/study/NCT05990166, August 2023.

Background: Conventional performance metrics do not indicate whether artificial intelligence (AI)-platforms and radiologists are right or wrong on the same examinations, nor whether their non-agreement is more often beneficial or harmful. Objective: To develop, illustrate, and apply directional metrics quantifying non-agreement between radiologists and AI-platforms against a common reference standard, and to adapt a fragility-based framework for quantifying the statistical robustness of these metrics. Methods: We defined Positive Non-Agreement (PNA, AI-platform correct, radiologist incorrect), Negative Non-Agreement (NNA, radiologist correct, AI-platform incorrect), and Net Non-Agreement (NENA, their difference). Metrics were illustrated in three hypothetical scenarios of 1,000 examinations each and applied to three previously published AI-radiologist studies on pulmonary embolism and mammography. Bootstrap confidence intervals were calculated, and a reclassification fragility index (FI) quantified the minimum case reclassification needed to remove statistical support for NENA's direction. Results: In illustration, NENA was positive in two scenarios and negative in one. In the three published studies, NENA ranged from -4.18% to +10.74% and decreased consistently with increasing radiologist experience (junior residents, +6.27%; consultants, +1.45%), a gradient driven primarily by declining PNA with modest increase in NNA. The reclassification FI was defined for all three illustrative scenarios (range, 0.10%-1.70%) and for six of the seven comparisons from these studies (range, 0.34%-4.31%), but not for the weakest comparison, whose confidence interval for NENA included zero. Conclusions: PNA, NNA, and NENA provide a directional, case-level complement to conventional performance metrics for AI-radiologist non-agreement. The accompanying reclassification FI offers a statistically grounded, generalizable method for quantifying their robustness.